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Abstract 

The  viscous  hypersonic  flow  past  the  leading  edge  of  a 
sharp  flat  plate,  whose  surface  is  parallel  to  an  oncoming 
uniform  flow,  is  analysed  on  the  basis  of  a  continuum  model 
consisting  of  the  Navier-Stokes  equations  and  the  velocity-slip 
and  temperature-jump  wall  boundary  conditions.  It  is  assumed 
that  the  model  fluid  is  a  perfect  gas  having  constant  specific 
heats,  a  constant  Prandtl  number,  <r,  whose  numerical  value 
is  order  unity,  and  a  normal  viscosity  coefficient  varying  as 
a  power,  w,  of  the  absolute  temperature.  Limiting  forms  of 
the  solutions  for  such  a  flow  are  studied  as:  (1)  the  free- 
stream  Mach  number,  M,  goes  to  infinity;  (2)  the  free-stream 
Reynolds  number  based  upon  the  distance  from  the  leading 
edge,  Rr  ,  goes  to  infinity;  and  (3)  the  ’Newtonian parameter', 
c=(Y-l)/(Y+l>.  where  y  is  the  ratio  of  the  specific  heats,  goes 
to  zero;  such  that  the  Newtonian  hypersonic  interaction 
parameter,  =(c5+wM2'2+u,)/Rl)1/2,  goes  to  infinity. 

Through  the  use  of  asymptotic  expansions  and  matching, 
it  is  shown  that,  for  «<1,  near  the  leading  edge, the  interaction 
of  the  thin  (thickness  ratio,  6,  going  to  zero),  viscous,  high 
temperature,  principal  layer  adjacent  to  the  plate  and  the 
external  flow  produces  three  distinct  strong-interaction 
(N=l/eM2  62<<l)(sub-)regimes,  for  which  the  quantity  6  is  of 
order  e^mM/tl+w)),  eVU+w^  and  eV2,  respectively,  as  the 
various  limits  are  approached. 

1.  Introduction 

According  to  ordinary  hypersonic  strong-interaction  theory 
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(OHSIT)  with  e  =  0(1),  for  flow  past  a  flat  plate  (Bush,[l]),  for 
«<1,  at  the  plate  surface,  there  is  a  thin,  viscous,  heat  con¬ 
ducting  layer,  v/hich  disturbs  the  external  uniform,  high 
speed  flow.  This  layer,  with  6,  the  layer's  thickness  ratio, 
much  less  than  0(1),  acts  as  an  effective  slender  'body',  pro¬ 
ducing  an  oblique  Rankine-Hugoniot  shock  wave  (SW),  whose 
distance  from  the  plate  is  also  of  0(6),  and  an  HSDT  inviscid 
shock  layer  (SL)  between  the  clearly  defined  outer  edge  of  the 
viscous  'body'  and  the  downstream  side  of  the  SW.  In  addition, 
it  is  found  that  there  must  be  a  viscous  HSDT-type  transition 
layer,  whose  thickness  ratio  is  much  less  than  0(6),  between 
the  viscous  'body'  and  the  inviscid  SL,  in  order  to  insure 
uniform  matching  between  these  two  layers.  This  theory  is 
shown  to  be  valid  for  )(=  M2+u/R jJ/2»1. 

The  Newtonian  hypersonic  strong-interaction  theory 
(NHSIT)  with  e«l,  is  formulated,  assuming  that  the  Newtonian 
equivalent  of  the  above  requirement  holds  true,  namely, 
X=e^WiM!t“/RLf  »  1. 

Clearly,  it  is  expected  that  the  NHSIT  should  produce  a 
regime  that  is  the  direct  equivalent  of  the  OHSIT  picture  just 
described.  However,  with  the  degree  of  freedom  that  the  add¬ 
itional  limit  of  e-^0  offers,  it  is  also  expected  that  the  NHSIT 
should  produce  additional  (sub-)regimes. 

In  the  following  sections,  it  is  shown  that,  indeed,  the 
NHSIT  does  produce  the  direct  equivalent  of  the  OHSIT  regime 
plus  two  additional  (sub- )regimes  that  owe  their  existence 
entirely  to  c's  approaching  zero. 

2.  The  Equations  of  Motion 

Consider  the  (two-dimensional)  flow  of  a  viscous,  com¬ 
pressible  gas  past  a  semi-infinite  flat  plate.  Let  Xj  =  Lx  and 
yj  =  Ly  represent  the  Cartesian  coordinates  parallel  and 
normal  to  the  plate,  respectively,  with  the  origin  of  this  coor¬ 
dinate  system  at  the  leading  edge  of  the  plate.  The  length  L 
is  to  be  chosen  so  that  x  is  of  0(1)  in  each  regime  that  is 
analysed.  The  velocity  components  in  the  Xj-  and  y, -directions 
are  u^UqqU  and  Vj=UqoV,  and  the  pressure,  temperature,  and 
density  are  P^PqoP,  T^TooT,  and  pj-poop,  where  u^,  p^,  T^, 
and  Pqq  are  the  velocity  in  the  Xj-direction,  pressure,  tem¬ 
perature,  and  density  in  the  undisturbed  region  upstream  of 
the  flat  plate. 

The  gas  is  assumed  to  be  a  perfect  one  (p=pT),  having  (i) 
constant  specific  heats,  cy  and  Cp  ,  with  y=(cp  /c  ),  such  that 
£=(y-U/(v+i)«i;  (  ii)  a  constant  Fr^ndtl  number  of  *order  unity 
(o-=const.  =0(1));  and  (iii)  a  normal  viscosity  coefficient  pro¬ 
portional  to  a  power,  to,  of  the  absolute  temperature  (p1=p00p=p00’Iw, 


with  ]/Z  ^  w<  1)  and  a  bulk  viscosity  coefficient  taken  to  be  zero. 

The  two-dimensional  Navier-Stokes  equations  including  the 
Fourier  heat  conduction  law,  assumed  to  be  the  governing 
equations  of  motion  for  the  flow  of  such  a  gas,  are 


where  M2=  (PqqU^)  /YRx>)*  the  square  of  the  free -stream  Mach 
number,  and  HltIpoo^xjL/jaoo)*  Reynolds  number.  The 
analysis  presented  here  is  for  M2»l  and  Rl»  1. 

3.  The  Principal  Layer 

According  to  existing  hypersonic  strong-interaction  theory 
for  flow  past  a  flat  plate,  at  the  surface  there  is  a  thin, 
viscous,  heat  conducting  layer,  which  disturbs  the  external 
flow.  This  layer,  designated  here  as  the  principal  layer  (PL), 
whose  outer  edge  is  given  by  y=6YJx),  with  6,  the  layer's 
thickness  ratio,  much  less  than  0(1),  has  as  its  proper  dis¬ 
torted  Cartesian  coordinates 

x  =  xb,  y  =  6yb  .  (3.01) 

The  expansions  for  the  flow  quantities  in  this  region,  carried 
out  in  the  above  distorted  coordinates,  are  taken  to  have  the  for 

u  =  ub+  ••• ,  v  =  6vb4  ••• ,  p  =(M262)pb+  ••• ,  T  =(eM2)Tb+  •••  ,p=(62/^)pb+**» 

(3.02) 

with  the  variables  i^  of  0(1). 

For  these  expansions,  in  which  it  is  implicit  that 


N  =  l/e  M2  62  «  1, 

the  leading  terms  in  the  equations  of  motion  are 


(3.03) 


where  X  =  e  1+tJM2c,,/RL&4  and  K^=  (e/62).  The  quantity  X=0(1), 
in  order  that  there  be  a  balance  between  the  inviscid  and  vis¬ 
cous  terms,  while  K,  must  be  >0(1)(K^«1  is  rules  out  as  not 
being  physically  realistic  for  this  flow  problem).  Note:  The 
combination  of  X=  0(1)  and  N«1  yields 

x  =£(34u)/2m2+W/r£2  »  1. 


From  Eq.(3. 04),  it  is  seen  that  there  are  two  different 
structures  for  this  layer  depending  upon  the  magnitude  of  the 
parameter  K^.  For  ^>>1,  the  y-momentum  equation  becomes 


8pb 

= 


and  Eq.(3.04)  corresponds  exactly  to  the  flat  plate  boundary 
layer  equations.  For  ^=0(1),  although  the  rest  of  the  equations 
are  of  flat  plate  boundary  layer  type,  the  y- momentum  equation 
is  not,  since,  in  this  limit,  the  convection  and  viscosity  terms, 
as  well  as  the  normal  pressure  gradient  term,  are  retained. 
The  PL  with  Kb»l  is  the  Newtonian  equivalent  of  the  OHSIT 
viscous  boundary  layer,  while  the  PL  with  K|?=0(l),is  a  special 
case,  for  which  there  is  no  OHSIT  parallel. 

To  complete  the  picture,  the  boundary  conditions  that  are  to 
be  applied  to  the  PL's  equations  of  motion  must  be  specified. 

The  boundary  conditions  at  the  wall  (y2  =  0),  taking  into  ac¬ 
count  velocity-slip,  thermal-creep,  and  temperature-jump,  are 


u,,o  =  C,{|i1(p1pl)'l/2f^}0  +|{p,(p1Tir1  , 

vl(0=-O»  Tl,o  ~  TI|W  =C2  {pi(PjPj)  ^  <Fy^,0'  O*05) 

where  Tj  w  is  the  temperature  of  the  plate,  and  0j  and  C2  are 
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constants  of  0(1),  which  depend  upon  the  momentum  and 
thermal  accommodation  coefficients  of  the  surface  and  the 
gas  properties.  Recasting  Eq.(3,05)  in  terms  of  the  PL 
variables  gives,  at  yb=0, 


Tv 

\C,  / b 


(2w-l)/2. 


V' 


}  “°»vu  =°» 


D  t(2W-1  )/2 

T  -T  -  XC*  f  b 
b,o  Ab,w  "  j^i/2  '  pb 


^}.o  =0‘ 


(3.06) 


with  T.  =  Tj  /c  M2Tqo.  Again,  two  different  structures  arise 
dependfng  on'the  parameter  K^.  For  K^»  1,  the  wallboundary 
conditions  reduce  to 


u,  =  v,  =  T  -  T  =  0  , 
b,o  b,o  b,o  b,w 

the  usual  non- slip  boundary  conditions.  On  the  other  hand, 
for  Kjj=0(l),  the  velocity-slip  and  temperature-jump  effects 
must  be  retained.  Even  for  K^=0(1),  the  effect  of  thermal- 
creep  is  negligible. 

Thus,  for  Kjj»l,  which  corresponds  to  6«c  1  ^2,  the 
equations  of  the  PL  take  the  form  of  the  flat  plate  boundary 
layer  equations  and  satisfy  the  non-slip  boundary  conditions 
at  the  wall;  while,  for  Kb=  0(l)(6=0(c  1/4)  ),  the  PL  equations 
are  of  flat  plate  boundary  layer  type  except  for  the  y-mo- 
mentum equation, and  satisfy  the  slip  boundary  conditions  at 
the  wall.  These  two  couplings  of  the  equations  and  wall 
boundary  conditions  for  the  PL  are  the  only  consistent  coup¬ 
lings  (cf.,  Pan  and  Pr obstein,  [2 ] ,  [ 3  ]  ;  Garvine,  [4], [5]  ). 

The  specification  of  the  boundary  conditions  at  the  outer 
edge  of  the  PL  cannot  be  made  explicitly  at  this  time,  since 
these  conditions  depend  upon  matching  with  the  solutions  for 
the  adjacent  layer. 

4,  The  Equations  of  Motion  for  the  External  Lavers 
The  analysis  of  the  flow  in  the  layers  external  to  the  PL  is 
most  conveniently  done,  assuming  the  PL  to  be  thicker  than 
any  of  these  external  layers,  in  terms  of  a  boundary  layer 
type  coordinate  system  oriented  with  respect  to  the  outer 
edge  of  the  PL. 

If  8j  =  Ls  and  nj  =  Ln  are  the  boundary  layer  type  coordi¬ 
nates,  with  the  Sj-axis  along  the  PL's  outer  edge  and  the  nj- 
axis  normal  to  this  outer  edge,  then  the  Cartesian  and 
boundary  layer  coordinate  systems  are  related  by 


X 

•J 


'■J 


x  =  \  cos  $(t)dt  -  n  sin$(s),  y=\  sin$(t)dt  +  n  cos  $(s),  (4.01) 


where  $(s)  is  the  angle  that  the  PL*s  outer  edge  makes  with 
the  free-stream  direction.  The  velocity  components  in  the 
Sj-and  nj-directions  are  qs  ^u^qg  and  qn|i=u00qn.  The 
velocity  components  of  these  two  coordinate  systems  are 
related  by 


u  =  q  cos$  -  q  sin$  .  v  =  q  sin$  +  q  cos$  . 
ns  nn  ns  nn 


(4.02) 


In  this  new  coordinate  system,  Eqs.( 2. 01) -( 2. 04) become 

,qs  ®qB  ®qs  §■  \  1  i 

°* p(irTr+\^‘T  Vw^Mih 


(4.03) 


where  h  =  1  -  $'n  . 


Under  the  assumptions  that  the  PL  has  a  thickness  ratio 
of  6  «  1,  and  the  thickness  ratio  of  any  external  layer  is 
<<0(6),  it  follows  that 


s  -*  x,  $  — ►  60  (sin$  -*•  60,  cos§-*l,  h-*l), 


(4.04a) 


where, 


since 


<f>-+  Y'h.  (4.04b) 

5.  The  Shock  Layer 

The  PL,  as  formulated  in  Sec. 3,  with  Kjj»l,  acting  as  an 
effective  slender  'body1,  produces  an  oblique  SW  and  a  SL 
between  the  'body*  and  the  SW.  The  structure  of  the  Newtonia 
SL  is  given  in  this  section, using  the  Newtonian  HSDT  (Cole, 
16])  as  a  guide  for  determining  the  orders  of  magnitude  of  the 
flow  quantities.  According  to  this  theory,  valid  only  for 
Kjj»l,  the  coordinates  of  the  layer  are 


s  =  s^,  n  =  (e6)nh  (5.01) 

and  the  flow  quantities  are  expressible  as 


qg  =  1  +  62uh  +  •••  ,  qn  =  M)vh  +  •••  , 


P  =  (M262)ph+ ...  .  T  =  (EM262)Th  +  -,  p=(l/e)ph+ -..,(5.02) 


with  the  variables  f^  of  0(1). 

For  this  layer,  then,  the  leading  terms  in  the  equations  of 
motion,  Eq.  (4.03),  are 


(5. °3) 

where  K,  =  (e/6*'+w). 

Since  \=0(1),  the  ratios  of  the  orders  of  magnitude  of  the 
leading  viscosity  and  heat  conduction  terms  to  thone  of  the 
inviscid  terms  in  Eq.(5.03)  are  Ofl/K^).  Thus,  within  the 
Newtonian  framework,  Eq.(5.03)  describes  either  an  inviscid 


shock  layer  (ISL)  or  a  viscous  shock  layer  (VSL),  depending 
on  whether  the  quantity  is  greater  than  or  equal  to  CXI). 

The  ISL,  with  K^»  1,  is  the  Newtonian  equivalent  of  the  OHSIT 
SL,  while  the  VSL,  with  Kh=0(l),  is  a  product  of  the  Newtonian 
approximation  that  e  -*■  0,  and  has  no  OHSIT  parallel.  (Note: 

W^^for  w<1» 

and  the  relation  K^>0(1)  does  not  violate  the  condition  1. ) 

To  complete  the  specification  of  the  SL,  it  is  necessary  to 
determine  the  shock  relations  that  are  to  be  applied  at  the 
outer  edge  of  the  layer,  n^  E=^h^sh^  *or  N«l,  anc*  show  that 
the  resulting  solutions  match  at  the  SL  inner  edge  to  the  PL 
solutions. 

6,  The  Shock  Structure  and  Shock  Relations 
An  analysis,  similar  to  the  one  performed  to  obtain  the 
shock  structure  for  viscous  hypersonic  blunt  body  theory 
(Bush,  [7]  ),  finds  that  either  PL-SL  combination,  with  N«l, 
and  Kb»  1,  supports  a  shock  structure  (SS)  which  is  made  up 
of  three  thin  layers:  (l)a  very  thin  exterior  layer,  in  whichthe 
order  of  magnitude  of  the  flow  quantities  is  characterized  by 
their  magnitude  in  the  free-stream;  (2)  a  relatively  thicker 
middle  layer,  in  which  there  is  dissipation;  and  (3)  a  thin  interior 
layer,  which  acts  as  a  transition  layer  between  the  middle  layer 
of  the  SS  and  the  SL. 

Since,  as  has  been  noted,  the  analysis  of  the  three  layers  of 
the  SS  for  NHSIT  follows  the  lines  of  that  given  in  [7]  ,  only  the 
proper  expansions  of  the  flow  quantities  and  the  limiting  forms 
of  the  equations  of  motion  for  these  layers  are  presented. 

6.1,  The  Uniform  Upstream  Region 

In  the  uniform  upstream  region  the  flow  quantities  are 

qg=cos$  =  l-62(<£2/2)+»**,qn=  -sin§  =  -604*”,  p=p=T=l  .  (6.01) 

b.c,  The  Exterior  Layer  of  the  Shock  Structure 
The  coordinates  of  the  exterior  layer  (EL)  are 


6=s  , 
e 


(6.02) 

and  the  leading  terms  inthe  expansions  for  the  flow  quantities  are 
q  =  cos§  +  62N1^°u  +  •••  =  {l-62(0z/2)+ •••  }+ 62N^cru  4  ••• , 

SO  6 
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q^=  -sin$  +  6N3^4  °v  +  *"  =  {-60  +  •••  }+  6N3^4  °ve  +  *** » 
p  =  1  +  N^4orp  4  ••• ,  p  =  p  4  •••  ,  T  =  T  +  ■"  . 


(6.03) 


The  first  integrals  of  the  equations  of  motion  for  this  layer, 
that  match  to  the  upstream  boundary  conditions  of  Eq.(6.01), 

are  f  ,0U  Z'3v  \“) 

V  <#>Pe=°.  Pc- T  --  0,  <iuU  [jevn  '  T  Fn  )J  =  °’ 

C  C  “ 

*■ ve+  X[l'Csr]=0’  «(T  -mx[iT“^e]=0,  (6.04) 


subject  to  the  restriction  of  N  >  0(e4,r),  This  restriction  must 
be  consistent  with  the  restriction  from  Sec. 5  of 
N  >  oo/m^-M/UM). 

6.3.  The  Middle  Layer  of  the  Shock  Structure 
The  coordinates  of  the  middle  layer  (ML,)  are 

‘=Bnt  — <6‘05) 

h  J  h  -* 

and  the  leading  terms  in  the  flow  quantity  expansions  are 
qe=cos§+62um4•••*l+62^n-^^+•,•,  q^3 -sin$  +  6vm+  •••  =6 


p=l+pm+...,p  =  (1/N)pm+-..,T=(1/N)Tm+-  .  (6.06) 

The  first  integrals  of  the  equations  of  motion,  determined  by 
matching  with  the  EL,  are 

(»pJ(vm^)^=°.Ptn-(1+pJTm=0.  0VX[Tm  F^0- 

^m+{lTm&]=0-  ^V''m»44C5r]=°-  <6-07> 

L  m-1  L  m  J 

6.4.  The  Inner  Layer  of  the  Shock  Structure 

The  quantities  in  the  inner  layer  (IL)  have  the  following 
representations: 


(6.08) 


=s.,  n=e  6^.(8.)+ 

|w(s^)4  e  m  4  • 


=  e6 


Nh<6h>+ 


q  =  cos$  4  6 

8 


=  14  62{W-(02/2)}  4€62u.4  •••, 


qn=(c6)v.4  —  ,  p=(l^)p.4  ••• ,  p=(l/eN)p.4  •••  ,T=(l/N) 


Q(s.)4cT.4..i 
(6.09) 


The  first  integrals  of  the  equations  of  motion,  from  the  ML-IL 
matching,  are  + 

_ _  9v.  ■ 

pi*VNh*+  0=piv* +  0=O»  PrQPi=  °»  *|( 


Q' 


8u.  i 

<jj  i 


377 


r,  8v.t 


8T.-I 

Iq“‘ 

a w  8n. 

L  iJ 


4  <p  (Q-$2)  =-  0  . 


(6.10) 


6.5.  The  Shock  Relations 

jprom  matching  with  the  solutions  of  the  IL  (as  n.-*-- co), 
the  boundary  conditions  at  the  outer  edge  of  the  SL,  N^(s^), 
are  determined  to  be 


■  T?  9T. 


uh.E+  4^k^).e  =-^vy(^),E=*‘ 

Vh,E=Nh'(Th,E//<^)J»  fy.E^^h.E*  Ph,E=  •  t6,11* 


where  fh  E=*h(sh,  Nh(s  )  ). 

Beforb  discussing  the  inviscid  shock  layer  regime 
(Kh>>l)  and  the  viscous  shock  layer  regime  (K^  =  0(1)  ), 
it  must  be  re-emphasized  that  the  SS  analysis  presented 
above  is  valid  only  for  K^>>1.  Therefore,  forK|:>  =  0(l), 
the  whole  concept  of  a  SW,  as  well  as  that  of  a  SL, 
must  be  modified. 

7.  The  Inviscid  Shock  Layer  Regime:  6«e 
For  the  inviscid  shock  layer  regime  (ISLR ), K^»l(Kh»l), 
the  primary  layers  of  interest  are  the  PL  (defined  by  Eqs. 
(3.04)  and  (3.06),  with  Kb»l)i  and  the  ISL  (defined  by  Eqs. 
(5.03)  and  (6.11),  with  Kh»l). 

It  is  possible  to  write  the  equations  for  these  layers  and 
the  boundary  conditions  at  the  SW  and  at  the  plate  in  imilarity 
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form,  as  is  done  for  OHSIT  in  [l],  if  N«1  and  Yb(xb)~xb^ 
and  Tb  w=  const.  However,  as  in  the  OHSIT  case, the  solutions 
of  the  and  ISL  do  not  match  directly  to  each  other,  and  a 
transition  layer  (TL),  intermediate  to  the  PL  and  ISL,  is  re¬ 
quired  to  provide  uniform  matching  for  the  entire  ISLR. 

Since  the  discussion  of  such  PL-TL-ISL  similarity 
solutions  is  presented  in  detail  in  [  1 )  ,  in  this  paper,  it  is 
considered  sufficient  to  present  the  equations  of  motion  for 
the  Newtonian  TL,  showing  the  dependence  of  the  flow  quanti¬ 
ties  in  this  region  upon  Kb»  1. 

The  quantities  in  the  TL  have  the  following  representations: 

s  =  st,  n  =  (e  6)Atnt,  (7.01a) 

q6  =  l+62atut  +  “*,  qn  =  M)Ptvt+  , 

p  =  (M262)p.  4-,  T=(eM262)0  T  +-,p=(l/e)(l/0  )p  4-, 

1  (7.01b) 

where 

6=  <*t=K^2  +",»  l,At=Pt=Kh1/'i+“)«  l,for  Kh»  1,  {  7. 02) 

and  the  variables  f^  are  of  0(1).  (Note:  It  is  also  true  that 

(e  M262)0t=  (e  M2)(e  6)2^24u),«e  M2,62at=(e  6)2^2+w)«  1.  ) 

From  Eq.(7. 01),  then,  the  leading  terms  in  the  equations  of 
motion  for  this  TL  are 


8.  The  Viscous  Shock  Layer  Regime:  6  =  Q(e^  b>^) 

For  the  viscous  lihock  layer  regirne  (VSLR),  Ki=Q(l) 
(Kb»l),  the  primary  layers  of  interest  are  the  PL  (defined 
by  Eqs.(3.04)  and  (3,06),  with  Kb»l),  and  the  VSL  (defined  by 
Eqs.(5.03)  and  (6.11),  with  Kb=0(l)).  (Note:  For  K^=0(1),  the 
thickness  ratio  of  the  ML  of  the  SS  is  of  the  same  order  of 


magnitude  as  that  of  the  SL,  namely,  0(e6).  ) 

For  this  regime,  it  is  not  possible  to  write  the  equations 
and  boundary  conditions  for  these  layers  in  similarity  form, 
due  to  the  additional  terms  in  the  shock  relations.  However, 
it  is  still  possible  to  show  that, since  the  VSL  equations  now 
contain  the  viscous  terms  that  are  required  in  the  TL  equa¬ 
tions  for  the  TL's  matching  with  the  PL  (cf.,Eq.(7. 03)),  the 
VSL  matches  directly  with  the  PL,  and  no  intermediate  layer 
is  required.  / 

9.  The  Merged  Layer  Regime: ;  6  =  Q(e  1'*') 

Consider  that  K|3=  0(1),  so  that  6  =  C^e^) « 1.  When  this  is 
the  case,  the  analyses  presented  for  the  SL  and  SS  cannot  be 
valid (cf.,  Sec. 5  and  6).  However,  for  K^=  0(1),  Eq.(3. 04), 
giving  the  equations  of  motion  for  the  PL,  contains  all  the 
terms  that  are  present  in  the  SL  and  SS  equations,  as  well 
as  those  in  the  viscous  boundary  layer  equations.  Therefore, 
it  is  concluded  that,  for  K|;)=0(1),  a  merged  layer  regime 
(MLR),  upstream  of  the  regimes  of  Secs.  7  and  8,  exists,  in 
which  the  three  above-mentioned  regions  are  replaced  by  a 
single  PL,  which  takes  on  all  of  their  properties. 

Strictly  speaking,  since  T  and  p  are  of  0(cM2)»l  in  the 
PL  of  the  MLR,  this  PL  now  includes  only  what  were  formerly 
the  SL  and  the  ML  and  IL  of  the  SS,  and  must  be  complemented 
by  an  EL,  similar  to  the  one  formulated  in  Sec. 6. 2,  with  Z  =  0, 
in  which  T  and  p  are  of  0(1).  From  the  matching  with  thiseEL, 
the  boundary  conditions  at  the  outer  edge  of  this  PL,  Y^x^), 
are  found  to  be 


Ub,E=1,Vb,E=  °»  Pb.E^b’  V^.E^^.E'^VW,1: 

(9. 01) 


The  boundary  conditions  at  the  wall  are  those  given  in  Eq. 
(3.06)  for  ^=0(1).  Further,  as  x^-*  0,  the  flow  quantities  at 
the  wall  should  also  approach  their  free-stream  values, namely, 


u,  =  1,  v.  =0,  p.  a=K, 
b,o  o  b,o  o  b,o  o  T> 


*  Tb,oo  “pb,o  o  "  °^b,oo  _fb^°j9° ^ 


An  examination  of  Eqs.(3.04),  (3.06),  (9.01),  and  (9.02) 
shows  that,  for  this  regime,  it  is  not  possible  to  write  this 
set  of  equations  in  a  general  similarity  form,  however,  it  is 
possible  to  show  that  local  similarity  can  exist  near  x,  =  0. 

Consider  that,  as  x^-*-  0,  the  coordinates  and  flow  quanti¬ 
ties  may  be  expressed  as 


J  * 

xK=^b*  yb  =  Yb^b>nb,with  4  “*  (A  =  const.);  (9.03) 


v 1  *  C  ub<  V  ♦  ••• .  v  W  V +  *"  *  pbe  V  W  V+-  • 

,  (9.04) 


b  b 
J 

-  c  v\/*/ 


P  D  */ 


ij  H*(nbM  - .  Pb=  KbeJbT  H*(nbM 


with  J  ,  J  ,  J  ,  J  ,  J^,>0. 
y  u  v  p  T 

Take  the  values  of  these  exponents  to  be 
1/2  <  Ju=2Jv=J  =JT=l/tl4u)s2/3. 

2/3  «Jy=(2«41l/2(Hw)<  3/4  (since  ]/2  «u<l).  (9.05) 

(Note:  The  value  taken  for  Jy  in  Eq.(9.05)  indicates  that  the 
effective  'body'  becomes  progressively  blunter  as  the  leading 
edge  is  approached,  rather  than  becoming  straight.)  Substitu¬ 
tion  of  the  above  expressions  into  Eq.(3.04)  (with  Kl=  0(1)  ) 
yields 


|yj;, 

t 


dU*- 

Vbdvi 


4 

KbL 


.  dR* 

|J_R7‘- J  T|,  -s - — 

T  b  y  'bdqb 


■]*± 
J  A* 


dV. 


=0» 


dU* 

MT-J  7J - 

1  T  b  y  'b  dqbJ 
dV*  - 


_  x  f  d  /H*w _ ) 

„  d’lb^b  dnb/  '  ’ 

r  dVIl  1  A  A  /  dVU\ 

[JvVb  -  Jy  %  ]  -  7^  (  3  (»T  } 

+  ^Tb  (!^(Hb“[JTU:-Jy%5]) 
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MdU*  ^  /  dU*  q 


dH* 
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THb-Jy%d^ 


{  i 


A*2Ku  o’  dq 


sr(» 


% 

wdHSv. 

3%^’ 


(9.06) 


t/(i+w) 


neglecting  terms  proportional  to  £V',TU,/.  The  corresponding 
boundary  conditions  at  the  outer  ( =  1)( derivable  from  Eq. 
(9.01)  )  are 


14 


j*  . 

d,E 


u:  „=  v,^  =  r 


* 

'b,E 


♦ 

b,E 


=  H* 


b,E 


=  0 


(9.07) 


Within  this  framework,  the  boundary  conditions  at  the  plate 
(Tlb=0),  with  respect  to  tangential  velocity-slip  and  no 
normal  flow,  become 


1  - 


XC, 


A*K 


V* 


^H*(2C0-l)/2 


(9.08) 


As  far  as  the  temperature -jump  condition  at  the  plate  is 
concerned,  there  are  two  such  conditions  that  render  the 
equations  of  motion  and  all  their  boundary  conditions  locally 
self-similar  as  0.  The  first  is  that  of  a  'hot  wall'  with 

w  —►  T^  w=  const,  as  0.  For  this  'hot  wall1  case,  the 

temperature-jump  condition  becomes 


T*  4 
b,w 


\C; 


A*K 


ui  j  . 

}  =  0  . 
d%  ,0 


(9.09a) 


second.is  that  of  an  'adiabatic  wall,'  with  Tu 
iSt.Kb,/(l?U,)  =6b1/'1+“  Hg,w=|bV(^)  Hg,0basW|b^O.  Fi 
'adiabatic  wall'  case,  the  temperature -jump  condition  i 


The  second, 
(const 
the 


(9.09b) 


Thus,  Eqs,(9. 06)-(9.  09 )  represent  the  'starting'  ordinary 
differential  equations  and  their  boundary  conditions.  Their 
solutions  provide  the  'initial  conditions'  for  the  parabolic 
system  of  partial  differential  equations  that  defines  the  con¬ 
tinuum,  thin  layer  model  that  has  been  postulated. 
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Fig,  1.  Schematic  Diagram  of  the  Newtonian 
Hypersonic  Strong-Interaction  Regimes  for 
Flow  Past  a  Flat  Plate. 
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13  abstract 

The  viscous  hypersonic  flow  past  the  leading  edge  of  a  sharp  flat  plate,  whose 
surface  is  parallel  to  an  oncoming  uniform  flow,  is  analysed  on  the  basis  of  a 
continuum  model  consisting  of  the  Navier-Stokes  equations  and  the  velocity-slip  and 
temperature- jump  wall  boundary  conditions.  It  is  assumed  that  the  model  fluid  is  a 
perfect  gas  having  constant  specific  heats,  a  constant  Prandtl  number, A"  ,  whose 
numerical  value  is  order  unity,  and  a  normal  viscosity  coefficient  varying  as  a 
power, ,  of  the  absolute  temperature.  Limiting  forms  of  the  solutions  for  such  a 
flow  are  studied  as:  (I)  the  free-stream  Mach  number,  H,  goes  to  infinity;  (2)  the 
free-stream  Reynolds  number  based  upon  the  distance  from  the  leading  edge,  R.  ,  goes 
to  infinity;  and  (3)  the  'Newtonian  parameter',  €  *(f- l)/(*+l) ,  where  t  is  trie 
ratio  of  the  specific  heats,  apes *  to  zero;  such  that  the  Newtonian  hypersonic  inter¬ 
action  parameter  ,%€  -  (C^M2'2  '/Rl)  ,  goes  to  infinity. 

Through  the  use  of  asymptotic  expansions  and  matching,  it  is  shown  that,  for 
ft)<l,  near  the  leading  edge,  the  interaction  of  the  thin  (thickness  ratio,  I  , 
going  to  zero),  viscous,  high  temperature,  principal  layer  adjacent  to  the  plate 
and  the  external  flow  produces  three  distinct  strong- interact  ion  (N“ 1 /€  M^S2^c  l) 
(sub-)  regimes,  for  which  the  quantity  S  is  of  order  < m (m>l/ ( !-»-•>) ) ,  61/0+1*), 
and  respectively,  as  the  various  limits  are  approached. 
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